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5^ ■ Abstract 
ON 

We study the Bloch theorem which states absence of the spontaneous current in 
interacting electron systems. This theorem is shown to be still applicable to the 
system with the magnetic field induced by the electric current. Application to the 
spontaneous surface current is also examined in detail. Our result excludes the 
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possibility of the recently proposed cZ-wave superconductivity having the surface 
flow and finite total current. 
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o\ : 1 Introduction 



Although the symmetry of the order parameter in high-T c superconductivity has not been 
identified yet, we have some experimental evidences which indicate the c?-wave superconduct- 
ing state. 0, @, H, % |], H 0, f| Since this type of the symmetry is qualitatively different from 
the conventional s-wave superconductivity, we can expect some new phenomena in high-T c 
sup er conduct i vity. 

Recently, some researchers pointed out the possibility of the d-wave superconductivity with 
broken time-reversal symmetry, which results from the s + id-wave state near the surface. 0, 
[T0| , [Tlj In addition, it was discussed that this state has a finite electric surface current. [0, [13 1 



In this case, we should remark that this state has the following properties: 
(1) The current flows without any external field, i.e., it appears spontaneously. 



(2) The total current is also finite. flT2 



On the other hand, the Bloch theorem states absence of the stationary current in the ground 
state with no external field. [13, [TJJ [TB|, [T7|, [TB| Then the above surface current seems inconsis- 



tent with this theorem. 

Motivated by this discrepancy, we study the Bloch theorem and its application to this state. 
The Bloch theorem has been examined so far within the framework of the first quantization, 
which seems disadvantageous to superconductivity. Furthermore, the theorem in case with the 
magnetic field induced by the spontaneous current has not been examined in detail. For these 
reasons, we present the second quantized version of the Bloch theorem, which can directly 



Submitted to J. Phys. Soc. Jpn. 



1 



cover the superconductivity, and also show its extension to the system with the induced 
magnetic field. 

This paper is organized as follows: In §2, we examine the Bloch theorem neglecting the 
magnetic field induced by the current, because the spontaneous surface current is reported 



in this situation. [12, [L3[ Effects of the induced magnetic field are discussed in §3, which is 



followed by the summary in §4. 



2 Bloch Theorem 

In this section, we first of all rederive the Bloch theorem from general viewpoint. Possi- 
bility of the spontaneous surface current is then discussed on the basis of this theorem. 

2.1 Definitions 

We start with the case of the chargeless (neutral) fermion, in which the current does not 
induce magnetic field. We also assume the absence of the external magnetic or electric field. 
Then, the Hamiltonian is given by 

» 2 

Here p = while m and fi represent the mass of an electron and the chemical potential, 
respectively. ^(r) means the electron field operator with spin-a, and U(r — r') is the 
electron-electron interaction, in which the pairing interaction can be included. The lattice 
and impurity potentials are represented by V{r). (To describe the presence of the surface, 
we put the potential V{r) infinity outside the sample.) 
We also introduce the current density operator 

Hr) = ^ £[*t(r)#Mr) - (p*t ( P ))^( P )]. (2) 

At this stage, we comment on the reason why this paper uses the second quantization 
in contrast to the previous papers, 0, [l7j which prove the Bloch theorem in the language 
of the first quantization. When superconductivity is examined, appropriate infinitesimal 
symmetry breaking field is formally necessary in order to select a special state from various 
superconducting states. In this case, the corresponding fictitious field Hamiltonian (Hp-p) 
has a number non-conserving form, and hence it is difficult to write down within the first 
quantization. We thus use the second quantization in order to avoid this problem. Indeed, 
Hpp for the superconductivity is given by 

H FF = J drdr'[/i(r,r')*}(r)*{(r') + h.c.} } (3) 

where h(r, r') is the infinitesimal symmetry breaking field. In what follows, though H FF is 
not written explicitly in eq. (|1|), the symmetry breaking field is always put zero after the 
system size is set infinity. 
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2.2 Ground state 

Firstly, we examine the system at T = 0. We assume that the ground state has a 
finite total current J^, 



J^p = / dr J^{r) 



dr<V| J(r)m = W\ E / dr*+ |V>> 7^ 0, (4) 



where J^(r) denotes a local current. In the following discussion, we show that this assumption 
of finite J ^ leads to a contradiction. 
Let us consider a trial state |L4] 



|0) = exp{WP ■ E / dr^t(r)r^ a (r)}\iP), (5) 

and compare (</>|i?|</>) with the ground state energy, (ip\H\ip). Physically, 8P in eq. (^j) is 
regarded as a fluctuation applied to \ip). If |^>) is the true ground state, it must be stable for 
arbitrary 5P, i.e., (i)\H\ip) < (<j>\H\4>). 

It is easily found that each of the (chemical) potential and the interaction terms gives the 
same expectation value in both states. On the other hand, the kinetic term gives 

(01 E / dr*t(r)|^<Mr)|0 = M £ / dr*t(r) & ^ P) V (r)M- (6) 
Using eq. (g), we obtain 

(0|tf|0> = + ■ Jv, + ^5P 2 , (7) 

2m 

where iV represents the number of electrons, 

N = f dr¥ a (r)* a (r)\i>). (8) 

cr ^ 

Now take 5P opposite to in direction and small enough so that the third term in the right 
hand side in eq. (|7|) is negligible compared with the second one. Then we find that \(f>) has a 
lower energy than 

(<f>\H\<j>) = - \6P - J* I < (V> W>- (9) 

Though |0) is not always the eigenstate of if, the variational principle guarantees that 
cannot be the ground state. This shows that the ground state of the chargeless fermions 
cannot have finite total current. We thus obtain the second quantized version of the Bloch 
theorem, which can cover the superconducting system. Moreover, we find that the theorem 
holds irrespective to the detail of the potential. We will use this result when the surface 
current is studied in section 2.4. 

We can extend this theorem to the system with magnetic impurities as well. Consider 

#mag = E/ drE«*(r " Ri)Si ■ ^l(r)(T a y^ a/ (r), (10) 

a, a' i 

where Ui is the potential of an impurity at Ri with spin Si, and cr represents the Pauli matrix. 
Then because of 

[tf mag , exp{i,5P ■ E / dr*t( r ) r «r ff ( r )}] = , (11) 

we obtain (0|i/ mag |0) = (tp\H maiS \ij;); the Bloch theorem is still valid even if magnetic impu- 
rities are present. 
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2.3 Finite temperature 

In the next step, we examine the thermodynamic stability of the spontaneous current. 
Introduce an orthonormal complete set of eigenstates an d define the density matrix p^ 

which gives the lowest free energy at T: 

Pip = E \i>i) w i{' t f>i\- ( 12 ) 
i 

Here Wi is the statistical weight which satisfies J2i w % = 1 an d < Wi < 1. (In many cases, the 
weight Wi is set as e~ Ei ^ T /Tr[e~ H ^ T ], where Ei represents the energy of the eigenstate 
We again assume that p^ gives a finite total current, 

J+ = f drJ*(r) = J drTr[p^J(r)] = Tr[p^ J dr*+(r)^* ff (r)] ^ 0. (13) 

In the following, we show that p^ is always accompanied by another density matrix which 
gives a lower free energy, i.e., cannot describe the most stable state. 

We consider a trial density matrix created by the operation of the exponential operator 
in eq. @: 

/^^EI&W&I' ( 14 ) 

i 

where 

\<Pi) = exp{i5P ■ £ / drtt* (r)r* ff (r)}|^). (15) 

It is found that {| also forms an orthonormal complete set. 

Since p^ and p^ have the same statistical weight, the entropy is also the same for these 
states, 

SV = - Tr [pV°gP</>] = - ^WilogWi = -Tr^logp^,] ee Sy. (16) 

i 

On the other hand, the expectation values of the energy of p^ (E^ = Tr[p^H]) and p^ 
(E^ = Tr[p^H]) are relating to each other in the form 



r SP • t> SP 2 

e^ = ^^l^i0i> = E^I^)+E^(^IE/ d ^( r ')( — - + ^-)*a(r)|^> 

_ J 7YI ZftTt 

% % % <y 

= E^ + SP-J^ + ^SP 2 . (17) 
2m 

From eqs. (^) and (0), we can compare the free energy for p^ (ee F^) with that for p^ 

F^ = E^- TS^ = F iP + 5P-J i> + ^P 2 . (18) 

For a SP being small so that the term with 0(5P 2 ) in eq. (|18"D can be neglected and having an 
opposite direction to Jf, the free energy of the trial state, F^, is smaller than F^ for arbitrary 
statistical weights. (Even if we put p^ the Gibbs state which must give the lowest free energy, 
F^ > F^ holds.) In conclusion, the state with a finite total current is not thermodynamically 
stable. 

Before ending this section, we briefly note Bohm's study for a finite current at T ^ 0. WM In 
his proof, the off-diagonal matrix elements of H for are dropped by using the properties 
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that (i) they are 0(8P) and (ii) from the perturbative analysis, they are found to affect the 
energy at most 0(5P 2 ). However, if there is degeneracy in the diagonal elements and, for 
example, = ((pj\H\(pj) {i ^ j) is realized, (ij)-element may affect the energy in 

0(5P). Hence this off-diagonal element cannot be neglected within the perturbation theory. 
On the other hand, since our proof does not rely on the perturbation theory, we do not meet 
this problem. 



2.4 Application to the spontaneous surface current 

Recently, the spontaneous surface current has been discussed in unconventional super- 
conductivity. We examine this problem on the basis of the Bloch theorem. 

Let us consider a semi-infinite system as shown in Fig. 1: Presence of the surface is 
described by putting V(x < 0) = oo. The surface may have a roughness. At the end edges of 
the system, we take the open boundary condition, so that we can choose arbitrary small 5P. 
(The system size is finally put infinity.) 

Suppose the presence of the spontaneous surface current as shown in Fig. 2(a). Then, 
what does the Bloch theorem state about this situation? The answer is that there must be 
another counter current inside the system which cancels the surface one as a whole (Fig. 2(b)). 
Indeed, we can find an example of this situation in the domain wall problem in unconventional 
superconductivity. |L9], ^] In this example, the local surface currents exist near both the sides 
of the wall. They, however, flow toward opposite directions to each other, thereby satisfying 
the Bloch theorem. At this stage, we emphasize that the presence of the counter flow purely 
results from the Bloch theorem and is irrespective to the origin of the surface current. 

In the next step, we discuss some specific systems in which the spontaneous surface currents 
are proposed. 

(1) As noted in the introduction, the surface current which is also finite in total was proposed 
in a (i-wave superconductivity with broken time-reversal symmetry near the surface. fTJL 



I3| Since the detail of the state, e.g., whether the time-reversal symmetry is broken or 
not, does not affect the Bloch theorem, we immediately find that it is not stable at all 
temperatures. 

(2) Sigrist et. al. also reported the surface current in an unconventional superconductiv- 
ity with broken time-reversal symmetry. [jl9| , |20[| They showed that the surface flow is 
canceled by the screening current inside the system resulting from the Meissner effect. 
Apparently, their results do not contradict the Bloch theorem. (Exactly speaking, the 
Meissner effect needs the charge of electron, which is neglected in this section. In §3, 
the Bloch theorem is shown to be still valid in the presence of the charge.) However, we 
have to note that the Meissner effect is not intrinsic origin of the back flow. The counter 
current has to appear even for the chargeless electrons where the Meissner effect does 
not work; the role of the Meissner effect is merely determining the distribution of the 
local current so as to exclude the magnetic field from the superconductor. 

Before going to the next section, we comment on the boundary condition at the end edges 
of the system. One may think that there is not stationary current under the open boundary 
condition (isolated system); however, when the thermodynamic limit is taken with appropriate 
symmetry breaking field, we can safely reach the bulk system in which the stationary current 
can flow. Furthermore, it is expected that the detail of the boundary condition does not affect 
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the presence of the surface current. Consequently, although we cannot examine the cases with 
arbitrary boundary conditions, the statement that the surface current state is not stable is 
believed to describe the property of the real system. 



2.5 Possibility of the spontaneous circulating current 

We examine the possibility of the circulating flow. In a real system, a circuit as shown in 
Fig. 3(a) is necessary in order to maintain the current. In this case, the total current passing 
through the line "A" in Fig. 3(a) is zero, which looks escaping from the Bloch theorem. Then, 
is the spontaneous current possible when the circuit is taken into account? 



Recently, Vignale extended Bohm's proof to systems with toroidal geometry. fl21fl The outline 
of Bohm's proof is as follows: Consider a large but thin superconducting ring with the width 
L ^ R (Fig. 3(b)). The ground state is assumed to have a finite total angular momentum 
as a result of the circulating current. The trial state is 

|0) = exp{in]T / drt>l(r)9t> a (r)}W). (19) 
Because of the single- value condition of \<j>), n must be an integer. We find 

(</>\H\t/>) = (i/>\H\i/>) + (V>|£/dr*t(r)^* CT (rM 

a J mr z 

+ (VlE/d^W^Mr)!^), (20) 



where L z is the z-component of the angular momentum operator. When the total current is 
Jip, the angular momentum is estimated as (L z ) ~ rriRJ^p. Then the second and third terms 
in eq. (|20"D are evaluated as nJ^/R and n 2 N/(2mR 2 ), respectively. Now, let us examine the 
case with = O(R). Because of N = O(R), the third term in eq. (|20|) becomes negligible 
compared with the second one for large R. Thus, when we choose n so as to be nJ^ < 0, it 
is found that \<p) has a lower energy than the assumed ground state. Consequently, the state 
that has a finite angular momentum due to a circulating current with O(R) cannot be the 
ground state. 

In conclusion, even if the circuit is taken into account, the spontaneous current which is 
proportional to the size along the circuit is not realized in the bulk system. 



Bohm proved for a large ring that the ground state must be in zero total angular momentum. [14 



3 Extension of the Bloch Theorem to the Case with 
Induced Magnetic Field 

In this section, we examine effects of the charge. Since the real electron has the charge 
e, the electric current induces magnetic field. Then, one may expect that the spontaneous 
current can be stabilized by the induced magnetic field. This section is devoted to study this 
possibility. 
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3.1 Definitions 

We summarize the definitions and equations used in this section. The Hamiltonian and 
the current density operator are given by 

H = ]P J dr*t (,■)[[ # " ^ (r))2 - „] + V(r))9 a (r) 

+ lY,f drd ^Ur)*Ur'Mr-tS)Vj(tS)V (r) (21) 

cr,cr' 

and 

J(r) = T,H(r)p^(r) - (p*t (r))<Mr)] - £. A (r) £ *+(r)* ff (r). (22) 

Here A(r) represents the vector potential. (We treat the vector potential classically in this 
paper.) We also define the spin magnetization operator 

M s (r) = - 9 /iBE^(^)f-'^W. (23) 

where g and //g represent the g-value and the Bohr magneton, respectively. 

The vector potential A(r) originates from the local spin magnetization iVf g(r) = (Afs(r)) 
and the electric current «7(t*). The magnetic flux density B(r) and the magnetic field H(r) 
satisfy 

V xA(r) = B(r) = ^(r) + M a (r) (24) 

and the Maxwell equation 

yxH(r) = J(r). (25) 

In eq. (|24"D, /x denotes the magnetic permeability of the vacuum. Under the Coulomb gauge, 
y ■ A(r) = 0, we find, from eqs. Q2"4D and (P5|), 

AA(r) = -/A)J(r)-vxM s (r). (26) 

Once the wavefunction (T = 0) or the density matrix (T ^ 0) is specified, A(r) is determined 
by eq. ( [26] ) with appropriate boundary conditions. 



The energy of the electron system is obtained as the expectation value of eq. (|21f) . To 
obtain the total energy, we have to take the magnetic field energy also into account, which is 
given by 

E u = - f drH(r) ■ B(r) = - f drA(r) • J(r). (27) 



2J y ' y ' 2 
In the last part of eq. (^71), we have put the surface integration of H(r) x A(r) equal to zero 
by taking the surface much larger than the sample. (When the sample is isolated, the current 
must circulate in it, which looks like a magnetic dipole moment at a place far away from the 
sample. Then, since the corresponding vector potential and magnetic field behave as ~ r -2 
and r~ 3 , respectively (r: distance from the sample), the surface integration of their product 
disappears for an infinitely large surface. [22]) 



Exactly speaking, there is the Zeeman energy by the spin magnetization 

E z = - [ drM s (r) ■ H(r). (21 



However, the following analysis does not treat this effect. We will give a brief discussion 
about this point later. (Note that the Zeeman energy by the orbital magnetization due to 
the spontaneous current is included in eq. fl2l"D through the vector potential.) 



3.2 Zero temperature 

Let us start with the case at T = 0. As discussed in §2, we introduce which has the 
local current, J^(r), the total current, the induced magnetic field described by A$(r) 
and the local magnetization, M^{r). We also introduce \<p) defined by eq. (§) with the 
local current, J<f>(r) = J^(r) +6J(r), the vector potential, A^(r) = A^(r) +5A(r) and the 
local magnetization, Ms</>{ r )- Since the exponential factor in eq. (f|) does not change spins of 
electrons, the magnetization should be the same between the two states, \ip) and \(j>). Then, 
from eq. (£§) for the sets (A^(r), J^(r), %W) and (A^r), J^(r), M s ^(r)), one finds 

A <L4(r) = —fi 5J(r), (29) 

where the explicit form of 8J(r) is derived from the expectation value of eq. (|22|) to be 

(JJ(r) = -(^1 E *Ur)**(r)\mP ~ e5A(r)\. (30) 



The term with 0(5P 2 ) is dropped in eq. Q3Tj). Substituting eq. ( |3"Q) into eq. ( |2"9"D and putting 
<5A(r) = <5Pa(r), we find that the equation of a(r) is independent of SP; 5A(r) and 5J(r) 
are therefore estimated as 0(SP). Keeping it in mind, we henceforth drop terms with 0(5P 2 ), 
by putting 5P small enough so that such terms can be neglected compared with the ones 
with 0(5P). 

Difference between and \<f>) in the field energy up to 0(5P) is 

E M <t, = \ j drA ^( r ) ■ J 4>i. r ) 

= \ j drA^r) • J^(r) + J dr5A(r) ■ J^(r) = E M ^ + 5E M , (31) 

where we have put the surface integrations of 5A{r) x (y x A^(r)) and A^,{r) x (y x 5A{r)) 
equal to zero as noted before. Substitution of Jf(r) into 5 Em gives 

6E m = 0j\ fdr5A(r)--J2^Ur)[p-eA i ,(r)]^ a (r)\^). (32) 
J m a 

where we have used the constraint y ■ 5A{r) = y • [A^(r) — A^{r)} = 0. 

Calculation of the expectation value of the Hamiltonian is performed in the same manner 
as in §2. Noting that A^(r) (A^r)) has to be used in eq. (plf) in evaluating (ip\H\ip) 
((4>\H\4>)), we obtain 

(<f>\H\<l>) = (ij\H\*P) + SP-(4>\ fdr-J2*Ur)\p-eA^r)]y a (r)\i>)-5E u 

J m a 

= {^\H\^) + Up-J^-5E m , (33) 
Ji> = WE / dr*t(r)-|p - e^(r)]*„(r)|V»>. (34) 



where 



From eqs. ([H]) and ([33|), the relation between the total energy of \ip) (E^) and that of \<f>) 
is 

E+ = E M4> + (<f>\H\<t>) = E^ + -5P ■ JV (35) 



e 
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When 5P is put opposite to J ^ in direction and small enough so that dropped terms with 
0(5P 2 ) are less important in comparison with the second one in eq. fl35|), \<p) always has a 
lower energy than \if)); this means that any state with a finite total electric current is always 
accompanied by a state with a lower energy, and it cannot be the ground state. 

The present theory does not deny the possibility that the Zeeman energy given by eq. ( p8|) 
stabilizes the spontaneous current state with finite net current. When this kind of mechanism 
really works, however, the system has a finite local spin magnetization, which should be 
observed experimentally. 



3.3 Finite Temperature 

Discussion at finite temperature in the presence of the induced magnetic field is almost 
the same as that in section 2.3. The difference is that and are now given by 

/ E i , = J2 i w i [(ij i \H( y A^)\ij i ) +E M ^}, 

where H(A q p i u^) means that A^.^ir) is used in eq. (|2T|). Each term in eq. (|36| ) can be 
evaluated in the same way as in section 3.2, and we again reach eq. (|17D, where J ^ in the 
present case is given by 

J* =5>iM£ / dr^Ur)-[p-eA A (r)}^ArM)- (37) 

We thus arrive at the conclusion that the spontaneous current which is finite in total cannot 
be stable thermodynamically even if the magnetic field induced by the current is taken into 
account, except for the case that coupling between the spin magnetization and the magnetic 
field stabilizes this state. 



4 Summary 

In this paper, we have examined second quantized version of the Bloch theorem, which can 
explicitly cover the superconductivity, and applied it to the recently proposed spontaneous 
surface current. We have also extended the theorem to systems with (1) the magnetic field 
by the spontaneous current and (2) magnetic impurities. 

Now, we summarize the results in this paper. 

(1) Any state with finite spontaneous total current is not stable at all temperatures, even 
if the magnetic field induced by the current is taken into account. This statement is 
correct unless the Zeeman energy by the spin magnetization stabilizes this state. 

(2) The recently proposed spontaneous surface current, which is also finite in total, is not 
stable at all temperatures. (Note that the Bloch theorem does not exclude the possibility 
of the surface current itself.) If a state with surface current is realized, the theorem states 
the existence of another counter current inside the system, which cancels the surface 
one as a whole. 

We will explicitly show in a subsequent paper that the surface current state really has 



higher energy than the one with zero total current. [23 
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Figure Captions 

Fig. 1: Schematic picture of the semi-infinite system. V(r) = in the right hand side in the 
figure and V(r) = oo in the shaded region. In the figure, z-axis is perpendicular to the 
surface of the paper. 

Fig. 2: (a) Assumed ground state, which has a finite surface current. In the figure and 
also in the following figures, thick arrows represent the electric current, (b) A possible 
surface current state having a counter flow inside the system which cancels the former 
in total. 

Fig. 3: (a) Circuit of spontaneous current, (b) A large ring with the width L and the radius 
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